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Abstract 

An asymptotic solution of the thermodynamics of a BCS superconductor containing spin 1/2 
and 7/2 magnetic impurities, obtained recently in Ref. [ij is exploited to derive the expressions for 
critical magnetic field Hc{T) . The credibility of the resulting theoretical equations, which depend 
on the magnetic coupling constant g and impurity concentration c, is verified on the experimental 
data for the following superconducting alloys: LaCe, ThGd and SmRh4B4. Good quantitative 
agreement with experimental data is found for sufficiently small values of c. The discrepancies 
between theoretical and experimental values of Hc{T) for larger values of c in case of LaCe and 
ThGd are reduced by introducing the concept of the effective temperature T, which accounts for 
the Coulomb interactions between the electron gas and impurity ions. 

At low temperatures, the critical magnetic field is found to increase with decreasing temperature 
T. This enhancement of the critical magnetic field provides evidence of the Jaccarino-Peter effect, 
which was experimentally observed in the Kondo systems like LaCe, (Lai_^Ce^)Al2 and also in the 
pseudoternary compounds, including Sni-^Eu^^MogSg, Pbi-^^Euj^MogSg and Lai.2-xEux>Mo6S8. 

The effect of an external magnetic field Ti on a BCS superconductor perturbed by magnetic 
impurities was also studied. On this grounds, by analyzing the dependence of superconducting 
transition temperature Tc on Ti of (Lai_xCej:)Al2, we have shown, that for certain parameter 
values, external magnetic field compensates the destructive effect of magnetic impurities. 

PACS numbers: 74.20.-z, 74.25.Bt, 74.25.Dw, 74.70.Ad 
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I. INTRODUCTION 



The effect of superconductivity has been discovered over 100 years ago^. Since then, 
many extraordinary properties accompanying this phenomenon, including perfect diamag- 
netism, zero-resistance, magnetic levitation due to expulsion of the magnetic field from a 
superconductor, flux quantization, Josephson effect and vortex state have been explained. 

The promising future applications of superconductors mainly concern electric power 
applications^ and generation of high, uniform magnetic field, e.g. for magnetic resonance 
imaging purposes^. However, the wide-scale applications of superconductors are limited, 
because superconductivity has been proven to be very sensitive to the destructive effect of 
an external magnetic field and the current density. In order to extend the applications of 
superconductors it is strongly desirable to enhance the values of the critical magnetic field, 
critical currents and superconducting transition temperature. These topics were of high 
interest during past few decades^"—. 

The superconducting critical temperature, as shown by Bednorz and Miiller-iS and later by 
Hosono group^^, can be raised by adding new elements to the antiferromagnetically ordered 
host system. Superconductivity appears at small value of dopant concentration x. Subse- 
quently, the transition temperature Tc increases almost linearly with x, and after reaching a 
maximum at optimal doping level Xopt, decreases and finally falls to zero. According to this 
scenario, the existence of a parent magnetic correlations is viewed to be an essential feature 
of high-temerature superconductivity. 

On the other hand, Matthias et al.^ discovered that the superconducting transition tem- 
perature of lanthanum decreases, when small amount (1 at. %) of the rare-earth magnetic 
impurities are added. It was shown, that the depression of Tc increases with the impurity 
spin value and not, as it was expected, with the impurity magnetic moment. 

The experimental studies carried out by Matthias revealed another extraordinary prop- 
erty of superconducting alloys, namely the possible coexistence of magnetism and supercon- 
ductivity. The coexistence of superconductivity and long-range antiferromagnetic ordering 
of the rare earth R magnetic moments was later discovered in RMogSes (R = Gd, Tb and 
Er)^, RRh4B4 (R = Nd, Sm and Tm)^ and in RMogSg (R = Gd, Tb, Dy and Er)^. A sim- 
ilar overlap between superconductivity and ferromagnetism was observed in ErRh4B4^ and 
HoMogSg^. 
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This discovery was surprising, since magnetism and superconductivity had been beheved 
to be mutually exclusive, because the internal magnetic fields generated in magnetically 
ordered systems are much larger than the typical critical fields of superconductors. How- 
ever, as independently predicted by Fulde and Ferrell^^ and Larkin and Ovchinnikov^^ , the 
superconductor may overcome the pair-breaking effect of magnetic field by forming periodic 
regions of superconductivity separated by domains of aligned spins. In such FFLO-state the 
order parameter is spatially modulated along the field direction. There are strongly exper- 
imental suggestions for the occurrence of FFLO state in some heavy-fermion compounds, 
e.g. in CeCoIn^^. 

Paradoxically, under specific circumstances, the external magnetic field can even en- 
hance (instead of depress) the properties of superconductors, e.g. the upper critical field 
Hc2 of SmRh4B4 increases below the Neel temperature Tn = 0.87 K->^. The latter can 
be understood on the grounds of the Jaccarino-Peter effect^^, in which the external mag- 
netic field compensates the antiferromagnetic exchange interaction generated by the con- 
duction electrons antiferromagnetically coupled to the spins of localized magnetic moments. 
The Jaccarino-peter compensation effect has been observed in several systems, e.g. in 

Euo.75Sno.25Mo6S7.2Seo.S^. 

In this study, we shall not investigate the FFLO state, since we assume a spatially uniform 
superconducting order parameter. In such the case, only the Jaccarino-Peter effect can be 
considered as the possible source of the superconductivity enhancement. Before going into 
details, let us briefly summarize theoretical investigations of the superconductors containing 
magnetic impurities. 

Early theoretical investigations of the problem of magnetic superconducting alloys were 
founded on perturbation theory. Nakamura^- and Suhl et al.— explained this effect by treat- 
ing the s-d interaction Vs-d''^ as an additive term in the total Hamiltonian, which perturbs 
a BCS superconductor—. Balseiro et al.— studied a BCS superconductor perburbed by 
magnetic impurities interacting via a nearest neighbour Heisenberg potential. The resulting 
phase diagrams comply qualitatively with experiment. 

The well-known Abrikosov-Gor'kov theory^^ (AG) of dirty superconductors explains the 
strong decrease in due to magnetic impurities and also predicts gapless superconductivity, 
confirmed experimentally by Reif and Woolf^. Disagreement with this approach is observed 
in a number of Kondo superconductors, e.g. hsii-xCex^^'r^^ LaCe and LaGd5i and PbCe 
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and InCe^. 

The AG theory was generahzed to describe an increase of Hc2 below the Neel tempera- 
ture Tn- To this end, Ramakrishnan and Varmai^ extended AG theory to the case, when 
the magnetic ions are present in large concentration. They solved numerically Eliashberg 
equations including the effects of phonons, spin waves and elastic scattering in order to 
estimate the variations of the pairbreaking parameter. However, a detailed comparison of 
their results with experimental data for the upper critical field of SmRh4B4 have been not 
performed. 

Hamaker et al.^'' found a good fit of the expression for the upper critical magnetic field 
in Machida's theory for antiferromagnetic superconductors to their experimental data on 
SmRh4B4 and obtained good quantitative agreement. However, some of the adjustable 
parameter values does not comply with other experimental findings. 

The effect of magnetic impurities on superconductivity is still under debate. Kozorezov 
et al.— have shown, that trace concentrations of magnetic impurities may also result in 
significant changes in nonequilibrium properties of superconductors. A comprehensive review 
of recent developments in this field can be found in Ref. |39 



These yet unresolved issues, as well as some shortcomings of the models presented above, 
motivate the present work. We investigate the thermodynamics of the Hamiltonian 

H^-^'^^ = H^cs + V^^'\ (1) 

where 

H-Bcs = Ho + Vbcs, (2) 

and 

Ho = ^^knka, (3) 

kcr 

with = £k — fJ', riko- = al^^a^.^, is the free fermion kinetic energy operator and 

^Bcs = -|A|"^ ^G'kk'4+«-k-«-k'-ak'+, (4) 

kk' 

is the Cooper pairing potential, whereas 

M 

y = /Ar-i^a,^,„, (5) 



represents the reduced s-d interaction. |A| denotes the system's volume and G^k' t^gsI, 
symmetric, invariant under k — )■ — k or k' — )■ — k' and nonvanishing only in a thin band close 
to the Fermi surface, viz., 

Gkk' = Gox(k)x(k'), Go > 0, 
where denotes the characteristic function of the set 

V = {k : n — 6 < e-i^ < H + 6}, 



2m 

In equation (jS]), Sza denotes the spin operator of the a— th magnetic ion, whereas 

0-^ = 5]^ (^k+ - '^k-) , 

kgP 

is the spin operator of a conducting fermion. M is the number of magnetic impurities, 
the number of host atoms. 

We assume the perturbation implemented by the localized distinguishable magnetic im- 
purities to be a reduced long-range s-d interaction, which involves only the z-components 
of the impurity and fermion spin operators (Eq. ([5])). The reason for this simplification is 
that the thermodynamics of the resulting Hamiltonian = Hq + Vbcs + V^^'^^ admits a 

mean-field solution, the accuracy of which improves with decreasing impurity density. Fur- 
thermore, this solution is thermodynamically equivalent to the one obtained for H with a 
Heisenberg type reduced s-d interaction 
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replacing V (Ref. |40|, Sec. 6.2.5). (Similarly, the thermodynamics of a classical supercon- 
ductor can be explained in terms of a reduced BCS interaction, whereas a gauge-invariant 
theory of the Meissner effect requires a more general pairing potential.) The reduced form 
of Vh, obtained by rejecting the sum Yl ^k' in the s-d interaction Vs-d and restricting the 
remaining sum X^^k to a thin layer, is an approximation which resorts to the fact that 

k 

spin-exchange processes, and not momentum exchange processes, are primarily responsible 
for the Kondo effect caused by Vg-d^^ and to the low-temperature regime applied. 

In section |TT] the system's free energy F{H, (i) is determined within the Bogolyubov in- 
equality. The resulting approximating formula for the free energy is used in sections IIIII 



and IIVI to derive the expression for the critical magnetic field He- The credibility of the 
theoretical expression for He-, which depends on the magnetic coupling constant g and im- 
purity concentration c, is subsequently verified on the experimental data for LaCe, ThGd 
and SmRh4B4. Good quantitative agreement with experiment has been obtained. At low 
temperatures, the critical magnetic field is found to increase with decreasing temperature, 
similarly as in some antiferromagnetic superconductors. 

In section |V] the effect of an external magnetic field on superconducting alloys is studied. 
The resulting expression for the superconducting transition temperature TdTi) is exploited 
in section |VT] in order to explain the compensation effect. 

The theory presented improves earlier developments in this field. Apart from an explana- 
tion of the reentrant behavior of superconducting alloys^, it provides an explanation of the 
Jaccarino-Peter effect and good quantitative agreement with the experimentally measured 
critical magnetic field of several superconducting alloys LaCe, ThGd and SmRh4B4. 

II. UPPER BOUND TO THE FREE ENERGY IN TERMS OF THE BOGOLYUBOV 
METHOD 

The full Hamiltonian of the system 

= Ho + Vbcs + V^^'^ (7) 
can be expressed in the following form in terms of mean-field parameters u, rj: 

i7W = /,W(zy,r/) + 4^\ (8) 

where 

h^'^\u,v) = h + h^l + ImN{u' - (9) 
h = Hbcs + i^(^z, = -gM{u - 7]) (10) 

C = ^7^E^- + ^^"VE^- (11) 

a a 

= -N-^ [9{a. - 5..) -uN]\ lN-^J2i9^. - vNY. (12) 

a a 

The Bogolyubov inequality 

F{H, + H,)<F{H,) + {H^)^^, (13) 



with Hi = h'-^\u,r]), yields 

/3) < + (i^^^)>,(M, • (14) 

The parameters u and 1] will be now choosen so that they minimize the free energy 
F(/i(^^)(z/,r/),/3) = -/3-MnTrexp [-/3/i(^^)(z/, r^)] , viz., 

gF(/^W(z/,r^),/3) _^ dF{h^^'){u,v),P) _^ 
du ' 9?7 

The explicit form of these equations is 

\ ^ TrS'^exp [-Phimp] . . 

N^""'^'^ N Trexp[-/3/.„„,] ' ^''^ 

V = ^ {cr.h ■ (16) 
Using equations ffT^ . f lT5]) and f lTB]) one obtains 

(^?"">„„, = ((Sf>,n, - (Sj;,,, ) . (17) 

Z \ imp/ 

The inequality Tr(py4^) > (Tr(py4))^ , valid for any bounded self-adjoint operator A and 
density matrix p, shows that (^H^^-^'j ^^j^^^ < 0. Hence, from equations (fT4l) . (fT7|) one obtains 
the relevant upper bound to the free energy 

F{H^^'\ /3) < F(/i(^^)(i/, r/), /3). (18) 

According to Eq. (fTSil we ascertain that the thermodynamics of the original system, charac- 
terized by H^^'^\ is almost equivalent to that of h^^^\ provided r] and u are the minimizing 
solutions of the equations f|T5l) . f|T6|) . The consequences of disregarding a term 
are discussed in section [IVl 

The two equations (fTSl) and ( fT6l) can be reduced to a single one for z/. The only require- 
ment is g > 0. The general form of Eqs. (JT5l) . (|T6i) is 



i^ = /i(^-^) + /2(^), (19) 
V = fi{^-V)- (20) 

Let (? > 0, then /2 > 0. Furthermore, 



V = 1^-/2(1^), 
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(21) 



which yields the equation for v. 

^^ = /l(/2M) + /2M, (22) 

where according to Eqs. ( IT^ . (1201) : 

/i(z/) = (MiV/3)-i|-lnTrexp[-/3/i(z/,0)], (23) 

/2(i/) = (iV/3)-i|-lnTrexp[-/3/.S;ij. (24) 

III. MEAN-FIELD THEORY OF li 

The form of the Hamiltonian h, given by Eq. (ITOl) is analogous to 

i^BCs(H) = ^0 + ^BCS - /^B^(T„ (25) 

describing a system of electrons with attractive BCS interaction in the presence of an ex- 
ternal magnetic field "H (/ir denotes the Bohr magneton). The explicit form of the system's 
free energy F{h^'^'^\h',T]), (3) can be therefore derived by exploiting the Bogolubov-Valatin 
transformational^ and the method developed in Ref. 44 for H^csC^)- 

This has been done in Ref. ^ for spin 5 = 1/2 and S = 7 /2 magnetic impurities perturbing 
the BCS-superconductor. For further investigation, let us recall the final form of the free 
energy F^^^ and briefly recapitulate the results of Ref. [l. 

The free energy is given by the following equation 



F(^)= min|p^|A| / \l-A'E''fs{^,E,^ji'^) - ^-Hn[2cosh{PE) + 2cosh{gmA'^)] 
+ M'c-\uf!f^ - i(/f y) + + Eo(A = 0) + p^5^}, S = 1/2, 7/2, 



(26) 



where F^^^ is the free energy of impurity subsystem, given by Eqs. f l33l) and fl35|) . E 



y^AM-^, pF denoting the density of states at Fermi level, 

mpp 
27r/?,^ 

whereas Eq{A = 0) denotes the ground state energy of free fermions. Two last terms in 
Eq. (126!) are the contribution to the free energy density from one-fermion states, lying outside 
V. 
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The system's state is characterized, according to equation (126|) . by the minimizing solu- 
tion, {Am, z/m}, of the following set of equations for the gap A and a parameter z/, describing 
the impurity subsystem 

A = IgoP 1' |/3(/3,i?,e,/fV))d^ S = 1/2, 7/2, (27) 



eg sinh(/3^M/fV)) 
M cosh{/3gMf^^\iy)) + cosh(/3Ek) 



- = A(/f V)) + /f V) = ^ + # V), (28) 

where 



/3(A E, e, /.) = r (29) 



sinh(/3E) 
cosh(/3E) + cosh{g/3Mf^''>{iy)) 

The properties of a superconductor with magnetic impurities can be determined by solving 
this set of equations, which is supplemented by the following condition for the chemical 
potential /x: 

^Tr(nkaPo) = n, (30) 

kcr 

where nko- = oJ.^a]^^ is the fermion number operator and n denotes the average number of 
fermions in the system. It has been shown in Ref. [Ij, that this condition, takes the form: 



E[l-^/3(/3,i5k,a,/f)] =n. (31) 



Equation f l3T|l resembles the BCS equation for /i and the properties of /s are similar to those 
of /bcs = tanh(/3£'k/2), e.g. both functions are odd in ^k- The solution of equation ( 13T|) 
is therefore exactly the same as in BCS theory, viz., fi = sp. Thus, we assume that in the 
low-temperature scale the following relations hold: 

Equations fl27|) and fl28l) clearly possess the solution A = z/ = for all values of /3 > 0. 
At sufficiently large values of /3 one finds also other solutions, viz., {A 7^ 0, z/ = 0}, {A = 
0, z/ 7^ 0}, {A 7^ 0, z/ 7^ 0}. Accordingly, we distinguish the following phases: 

— paramagnetic phase P with {Am = 0, z^m = 0}, 



— unperturbed superconducting state SC with {A^ 7^ 0, z^m = 0}, 
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— ferromagnetic phase F without bound Cooper pairs and {Am = 0, z/m 7^ 0}, in which 
impurity spins tend to ahgn opposite to those of conduction fermions (cf. Eqs. f lTOj) 



— intermediate phase D in which superconductivity coexists with ferromagnetism and 



We define the following temperatures corresponding to the respective phase transitions 

— Tc, 2nd order transition SC — )■ P. 

— Tpp, Curie temperature of 2nd order transition F ^ P. 

— TscD, 1st order transition D — > SC. 

— TpD, 1st order transition D ^ F. 

— TscF, 1st order transition SC — F. 

The set of Eqs. (!27|) . ( 128|) were solved numerically for different values of g, M, 6 and CqPf 
and solutions, minimizing free energy are exploited in next section to determine the critical 
magnetic field. The parameters were adjusted to fit the experimental data for LaCe, ThGd 
and SmRh4B4. These results are discussed in section HVl 

The free energy (Eq. [22]), as well as the equations for A (Eq. and u (Eq. 

strongly depend on the value of the impurity spin S. In the present work we study the effect 
on superconductivity of the following magnetic ions: Ce {S = 1/2), Gd and Sm (5* = 7/2). 
For S" = 1/2 one obtains 



and (HI])). 



{Am7^0,z/m^0}. 



M 



imp 



/3-Mn Tr exp [-/^CJ = -Z?"' E ^r exp [-^hH] 



(33) 



a=l 




with 




cg_ 
M 



ta.nh.{(3gi'). 



(34) 



Accordingly, for spin 7/2 impurities 



imp 



-M/3-Mn2 exp[-24:g^ f3N-^] cosh{7 f3gu) +exp[-12g^f3N-^]cosh{5f3gu) 
+ exp[-4g'^PN^^] cosh(3/35(z/) + cosh{Pgu) + ^cg"^, 



(35) 
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where 




MR 



7exp[-24^^/3iV"^] sinh(7/3c/z/) + 5 exp[-12^2/3iV"^] sinh(5/3^z/) 



(36) 



+ 3exp[-4^2/3A^-^] sinh(3/3^z/) + sinh(/3^z/) 



and 



R = exp[-24c/2/3iV"^] cosh{7(3gu) + exp[-12c/2/3A^-i] cosh{5/3gu) 
+ exp [-45(2/3 A^"^] cosh{3l3gu) + cosh{f3gu). 



The complexity of the free energy F^^^ and functions /i, /2, /s increases with the impurity 
spin value. It follows, that the impurity spin is the key factor affecting the thermodynamics 
of superconducting magnetic alloys. This conclusion is complementary with the fundamental 
experimental observation made by Matthias et al.-. 

IV. THE CRITICAL MAGNETIC FIELD 

The critical magnetic field Hc<s> forcing a system to undergo the phase transition from the 
$ phase to paramagnetic (normal) phase P is given by the equation 



where po denotes the vacuum permeability, Fp and F$ denote the free energy of the P and 
$ phase respectively. 

The free energy of the normal state can be obtained from equation (1261) with A = and 
z/ = 0, which yields 

— for spin 1/2 impurities 



H^^ = {Fp - F$)/2/io 



(37) 




Hn2 + -cg^ 
2 ^ 



(38) 



+ Eo(A = 0) + pp5^, 



for spin 7/2 impurities 



F^}^ = -2pp\A 





(39) 



+ Eo{A = 0) + pp5\ 



11 



The critical magnetic field required to suppress superconductiviy (phase transition SC — )■ P) 
is given by the equation f l37|) with Fsc replacing It will be denoted as usual by He. The 
expression for Fsc results from equation with {A 7^ 0,z/ = 0}: 



SC 



-Pf\M 



2 ^ 



-2/3~Mn2cosh(/3E) d^ 



(40) 



+ -c/ + Eo(A = 0) + p^5^ for S = 1/2, 



and 
F^sc =Pf\M 



-5 



A"^- Vj^^^^(/3, - 2/3-^ In 2 cosh(/3E) 



- M/?-^ In 2 exp[-2A/3cg^M-^] - exp[-12 ^cg'^M'^] - exp[-4/3cg^M-^] + 1 
+ Icg^ + Eo{A = 0) + pp6\ for ^ = 7/2. 



(41) 



The credibility of the above theoretical expressions will be now verified on the experti- 
mental data of LaCe, ThGd and SmRh4B4. The parameters g, GqPf, 6, M were adjusted to 
provide the best possible fit of the theoretical critical magnetic curves to the experimental 
data. In order to perform a fit, the set of equations ( I27|) . (!28|) was solved numerically for 
fixed values of the parameters g, Pf|A|, GqPf, S, M. The minimizing solution of this set 
of equations was subsequently substituted to the expressions for He or Hcd depending on 
which of the system's states {SC or D) possess the smaller values of the free energy, i.e. 
which of them is thermodynamically stable. 

The critical magnetic field of LaCe and ThGd are depicted in Figs. [1^, [2^ and [31 The 
agreement with experiment is satisfactory. The discrepancies between theoretical and exper- 
imental data increases with decreasing temperature and increasing impurity concentration 
in case of LaCe and ThGd. However, these discrepancies can be reduced by introducing the 
concept of the effective inverse temperature /3 = 1/ ksT, which is related to the system's real 
temperature by the following formulae 

/3(/3,7)=7"'tanh(/37). (42) 

The effective temperature results by averaging (over the impurity positions) the single parti- 
cle equilibrium density matrix of a quantum particle in a field of randomly positioned wells, 
representing the screened Coulomb potential at each impurity site^. Furthermore, it has 
been shown, that the 7 is of the form 

h 



7 



'■^/Mu 



21^1 
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with U2 denoting the 2nd derivative at well's minimum. In the present work 7 will be treated 
as the adjustable parameter. 

Theoretical curves of the critical magnetic field of LaCe and ThGd superconducting alloys 
with /3 replacing /3 in Eq. ( 137|) are given in Figs.[T|D,[2]D. It is clear, that the application of the 
effective temperature improves the agreement of the given model with experimental data for 
LaCe and ThGd, proving that the Coulomb interactions apart from Coulomb interactions 
are also important factor of the superconducting alloys theory. The effective temperature has 
been also proven to be the crucial in description of the doping dependence of superconducting 
transition temperature Tc in high-Tc superconductors^. 

In case of SmRh4B4, the concentration of magnetic impurities (Sm'^"'" ions), c ~ 11% is 
much larger than in case of LaCe and ThGd. Accordingly, the agreement with experimental 
data for SmRh4B4 is only satisfactorly and can not be improved by the application of the 
effective temperature. This suggests, that the term {^//^^^^^^^^j can not be disregarded for 
sufficiently large c and shall be included in the free energy F(h^^'^\(3) computation or the 
dependence of the Ginzburg-Landau parameter k on T shall be taken into account. In the 
case of large c, the magnitude of exchange interaction between conduction fermions and 
magnetic ions may presumably exceeds the magnitude of the Coulomb attraction between 
magnetic ions and conduction fermions. 

The solutions of the set of equations for A and u under varying temperature and impurity 
concentration of Ce and Gd are depicted in Figs. HHH 

The graphs, depicted in Figs. [1^ and |2|l show an increase of the critical magnetic field at 
very low temperature scale and for sufficiently large values of impurity concentrations, e.g. 
for c = 2.0% Ce at Fig. [TJi and for c = 0.20% Gd at Fig. E]). 

This enhancement of superconductivity also increases with impurity spin and impurity 
concentration, since at extremely low temperatures, the following formula 

H^d{c = 0.20%Gd) > H^d{c = 0.10%Gd) (43) 

is valid. 

The above observation results from the following fact. If the solution {A 0,u ^ 0} 
minimizes the free energy, then the values of Hcd are larger than H^. Accordingly, one 
obtains an increase of the critical magnetic field below Tscd, since at this temperature, 
the system undergoes a phase transition SC — t- D and Hcd becomes equal to He- This 
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conclusion may be recognized as incompatible with physical intuition, suggesting that the 
external magnetic field should gain the perturbative effect of magnetic impurities. As a 
result HcD should possess smaller values than H^. 

An increase of the upper critical magnetic field Hc2 = k\/2Hc with k = A/^, denoting 
the Ginzburg-Landau parameter, has been observed in the following materials: SmRh4B4, 
GdMogSg, TbMoeSg, Sni_j,.Euj,.Mo6S8, Pbi„2,-Eu2:Mo6S8, Lai.2-xEu^Mo6S8'^'^. These experi- 
ments confirm the validity of inequality (H3!) . 

Fischer et al.— has pointed out that the superconductivity enhancement, represented by 
an increase of the critical magnetic field is a result of the Jaccarino-Peter effect^. This 
phenomenon may occur in the Il-type superconductor, in which the magnetic moments of 
the impurities are antiferromagnetically coupled to that of the conduction fermions. This 
interaction generates an exchange field Hj, which acts on the spins of conduction electrons 
equivalently to an applied magnetic field, viz., breaks the Copper pairs. However, the 
negative sign of the coupling between the magnetic moments and the conduction fermions 
spins, determines the direction of Tij to be opposite to that of Ti. Thus, an applied magnetic 
field will be compensated by an exchange field, since the net magnetic field "Ht is given by 
Ti — l^jj. A given compound displays superconducting properties as long as the following 
relation holds 

- Hp < < Hp, (44) 

where 



Hp = J^^- (45) 

V Xp - Xsc 

Xp and xsc denotes the magnetic susceptibility of the normal and superconducting state. 
Hp, defined by Eq. f HS]) is the Chandrasekhar-Clogston limiting paramagnetic field4^i^. 



The Jaccarino-Peter effect has been observed experimentally in Euo.75Sno.25M06S7.2Seo. 
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during the investigation of the upper critical magnetic field. In particular, at low temperature 
scale, three subsequent phase transitions has been observed with increasing value of an 
external magnetic field, i.e. SC — )■ P — > SC — > P. 

In conclusion, it is worth to point out, that the magnetic impurities were proven to 
limit the superconductivity, but on the other hand, under some specific conditions, they 
help the superconducting system to overcome the destructive effect of an external magnetic 
field. Furthermore, the interplay between superconductivity and magnetism is believed to 
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FIG. 1. The critical magnetic field of LaCe under varying Ce concetration for real temperature 
T and parameter values given in Table [J (a) and for the effective temperature T and parameter 



values given in Table lU (b). The points are experimental data from Ref. 
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be a possible mechanism of high-Tc superconductivity^^, since the undoped state of cuprate 
superconductors is a strongly insulating antiferromagnet. The existence of such a parent 
correlated insulator is viewed to be an essential feature of high temperature superconduc- 
tivity. 

In the above discussion we dealt only with the critical magnetic field. To fully judge, 
if the Jaccarino-Peter compensation may occur in the superconducting alloys, one should 
study the effect of external magnetic field 1-L on such system. This will be done in the next 
two section, in which we first study the effect of "H on the superconducting alloys and then 
apply the resulting formulas to investigate the dependence of the superconducting critical 
temperature Tc on "H. 



V. THE EFFECT OF AN EXTERNAL MAGNETIC FIELD ON THE SUPER- 
CONDUCTING ALLOYS 

In the presence of an external magnetic field the system's Hamiltonian (Eq. ([7])) should 
be supplemented by the additional terms, describing the interaction with a magnetic field. 
This yields^: 

{%) = H^''^ - ^g'liB-H Sza - f^BUa,, (46) 

a 

where g' is the modified Lande factor—, /iB denotes the Bohr magneton, Ti = goH is the 
effective magnetic field at each impurity site. 
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FIG. 2. The critical magnetic field of ThGd under varying Gd concetration for: real temperature 
T and parameter values given in Table [J (a) the effective temperature T and parameter values 



given in Table HIl (b). The points are experimental data from Ref. 



52. 




o!o 0^5 l!o 2,0 2.5 Zfl 



FIG. 3. The upper critical magnetic field Hc2 of SmR 
Table HI The points are experimental data from Ref. 
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14B4 for the parameter values given in 
Hc2 is related to the thermodynamic 



critical magnetic field He by the expression, Hc2 = k^/2Hc. The Ginzburg-Landau parameter k 
has been treated as the additional adjustable parameter used to fit experimental data and assumed 
to be independent in T. The value of k providing best fit is k = 1.51. 



The additional electron cr^ and impurity spin Sz operators in Eq. ( I46l) are separated. 
This allows to use thermodynamic equivalence of Hamiltonians and h^^'^^u,!]), proven 
in Sec. [TTl to describe the thermodynamic properties of a BCS superconductor containing 
the magnetic impurities and in the presence of magnetic field in terms of the Hamiltonian 
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1 2 3 4 5 

T[K] 



1 2 3 4 5 

T[K] 



FIG. 4. The temperature dependence of the gap parameter A(T) for LaCe and the solution of 
Eqs. dlZl), dig) with 5=1/2, in which {A / 0, = 0} (a) and {A / 0, / 0} (b). The parameter 
values are given in Table HI 




T [K] T [K] 



FIG. 5. The temperature dependence of the gap parameter A(T) for ThGd and the solution of 
Eqs. (IlTD, (USD with S = 7/2, in which {A / 0, / 0} (a) and {A / 0, = 0} (b). The parameter 
values are given in Table HI 

h^^\l-L)^ which has the following form: 

/,(M)(^) = ^) _ J2 S,^ - fi^-Ha,. (47) 

a 

Accordingly, the electrons and impurities are described by the Hamiltonians h and himp 
respectively, 

h = Hbcs + (/« - ftB'H)az 
= Hbcs - {gMjf^ + fi^n) a, (48) 
= -ffBCs — Aa^, 
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FIG. 6. The temperature dependence of the parameter i^(T) for the solution of Eqs. (j27p . ()28p . in 
which {A 7^ 0, 7^ 0} for LaCe (a) and ThGd (b). The parameter values are given in Table HI 



a 




T [K] T [K] 



FIG. 7. The gap parameter A(T) at the effective temperature T and under varying impurity 
concentration for LaCe (a) and ThGd (b). The parameter values are given in Table HIl 



-(2 
za 



" " (49) 

= 9CY.S^. + \g'N-'Y.Sl^- 

a a 

The form of equations fj47|) and fH9l) is very similar to Eqs. ffTOl) and ffTTl) . It follows 
that, in order to include the effect of an external magnetic field on the free energy of a 
BCS superconductor containing magnetic impurities it suffices to perform the following 
substitutions in Eq. (1261) : 

K — 7- —A = K — i^bH, (50) 
y^C = ^- ^a'f^B-H. (51) 
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Alloy 


c [%] 


M 


9 [^] 


5 [eV] 


GqPp 


LaCe 


0.00 





0.000 


0.01 


0.3030 


LaCe 


0.55 


2 


0.080 


0.01 


0.2890 


LaCe 


1.00 


4 


0.095 


0.01 


0.2800 


LaCe 


1.50 


5 


0.100 


0.01 


0.2650 


LaCe 


2.00 


6 


0.105 


0.01 


0.2500 


ThCd 


0.00 





0.00 


0.01 


0.2200 


ThGd 


0.10 


3 


0.11 


0.01 


0.2090 


ThGd 


0.20 


5 


0.18 


0.01 


0.1910 


SmRh4B4 


11.11 


7 


0.0625 


0.0106 


0.2540 



TABLE I. The parameter values exploited during the numerical analysis of Eqs. ()27p . ()28p and 

dSZl). 



Alloy 


c[%] 


M 


9 [VeV] 


5 [eV] 


GqPp 


7 [ 10"^ eV] 


LaCe 


0.00 





0.000 


0.01 


0.3030 


0.00 


LaCe 


0.55 


2 


0.080 


0.01 


0.2900 


1.20 


LaCe 


1.00 


4 


0.095 


0.01 


0.2810 


1.30 


LaCe 


1.50 


5 


0.100 


0.01 


0.2655 


1.50 


LaCe 


2.00 


6 


0.105 


0.01 


0.2520 


1.50 


ThGd 


0.00 





0.00 


0.01 


0.2200 


0.00 


ThGd 


0.10 


3 


0.11 


0.01 


0.2090 


0.20 


ThGd 


0.20 


5 


0.18 


0.01 


0.1915 


0.30 



TABLE XL The parameter values exploited during the numerical analysis of Eqs. (I27p . ()28p and 

m- 

Accordingly, the set of equations for the parameters A (Eq. (1270 ) and v (Eq. (!28l) ) in the 
presence of the external magnetic field take the form: 

A = ^Gop^/' §/3(/3,i?,e,/fHC),H)de, (52) 

V = h{P,E,ij!f\C),H)) + if\C,H), S=l/2, 7/2, (53) 
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where 



/i(/3,E,e,/f^(C),H) 



eg 



sinh 



/3(^?M/f)(C) + /iBH 



M 



cosh 



/3(/3,i?,e,/f^(C),^) 



(3(gMf^'\C)+fiBn) +cosh{l3E) 
sinh(/3E) 



(54) 
(55) 



cosh(/3E) + cosh ^(3 (^gMf^^\() + fi^n ^ 

Functions /2 (C) are given by equations ( |34l) . (!36l) . with ^ replacing u. 

The free energy of the BCS superconductor perturbed by magnetic impurities and in the 
presence of an external magnetic field H then reads 



+ 2 cosh 



-5 



/3(gMff\0+f,Bn) ] 



2„-l 



d^ + M'c 



+ ^S(^) + ^o(A = 0) + pp6'}, S = 1/2, 7/2, 



where SH) are given by Eqs. fl5^ and fl5Sl) after substitution C, —)■ u. 



VI. CRITICAL TEMPERATURE 

The phase diagrams of a BCS superconductor perturbed by magnetic impurities depicted 
by us in Ref.|l|show that, the phase transition from the normal (Non SC) to a supercondcting 
state can be of the first or second order, depending on the value of the magnetic coupling 
constant g. The next two subsections are concerned with computation of the transition 
temperature Tc{l-L) for first and second order phase transitions. 



A. Second order phase transitions 

According to section llllt Eq. (127|) for the solution {A 7^ 0, z/ = 0} reduces to the BCS 
gap equation 

Abcs = ^G'oPf/ ^tanhf^/3EBcsV^' (^7) 
Ebcs = \/e + A|cs- 

The transition temperature Tc^'"^'^ in BCS theory, is defined as the boundary of the region 
beyond which there is no real, positive Abcs satisfying Eq. (I57j) . Below Tc^'^^^ the solution 
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-^Bcs 7^ minimizes the free energy and the system is in superconducting phase. Therefore, 
Tp^^') can be obtained from Eq. (1571) with Abcs = 0, which yields^°: 



(58) 



It should be possible to estimate the change in Tc^'^^\ since the density of states enters 
exponentially in Eq. fISH]) . However, significant deviations from Eq. f l58l) were observed 
experimentally for a number of superconductors containing magnetic impurities. This inad- 
equacy of Eq. (!58l) is most distinct for large values of impurity concentration. BCS theory 
is therefore incapable to describe the superconducting alloys. 

Expression for transition temperature of a superconducting alloy in the presence of an 
external magnetic field for 2nd order phase transition can be computed analogously as in 
BCS theory. To this end, it suffices to put A = in Eqs. (!52l) and (l53l) . Thus, one obtains 
the following set of equations for Tc = l/{k(3c): 

GqPp 



sinh(/3,|e|) 



1^1 cosh(/3,|e|) + cosh /3e (gMfPiQ + 



eg 


sinh 


[gMf^'\Q + 






^cosh 


/3,(^?M/f)(Cc)+/iB^) 


+ COsh(/3e|e|) 



(59) 



(60) 



where = z/^ - \g ^g'pB'H, J^c = ^{Pc), S = 1/2, 7/2. 

Numerical analysis shows that in the low-temperature scale I'dT) is almost independent 
in T, viz. I'dT) ^ z/(0) = cg/M. Accordingly, the set of equations ( l59l) . ( l60l) is solved under 
the assumption, that = cg/M. 

The resulting solution for Tc('H) for S = 1/2, 7/2 under varying impurity concentration 
is depicted in Fig. |8l The solution for Tc('H) for small c is similar to the numerical result 
obtained by Sarma^ of the system described by the Hamiltonian Hs = Hbcs + /^b^ctz- His 
result for T^^^^\n) agrees qualitatively with Tc{'H) graphs depicted in Fig. [HI since the 
expression for Tc^'"^\'H) obtained in Ref. 55I is of the similar form to Eq. f l5I?]) with = 0. 

The Tc{'H) graphs depicted in Fig. [H^ show two phase transitions for sufficiently small 
impurity concentration, i.e. c G (0 at. %, 0.10 at. %) and for l-i G (5.65 T, 7.8 T)). The first 
phase transition (P — t- SC) occurs at Tci('H) and the second phase transition (back to normal 
state) appears at Tc2{'H)- In case of LaCe, only one phase transition is present {SC — t- P) for 
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higher concentrations, but superconductor overcomes larger values of an external magnetic 
field. Furthermore, Tc{'H) initially increases with "H. 

The form of denominator on the right hand side of Eq. (!59|) suggests that, the perturbative 
effect of magnetic impurities can be compensated by an external magnetic field. It follows 
from the fact that /2 (Cc) is odd function in Cc and from a definition of C,c parameter, 
which approaches negative values for sufficiently large "H. Thus, the magnetic field intensity, 
required for the full compensation of the perturbative effect of magnetic imurities on a BCS 
supercondutor has the form: 

W. = -^/f'(Cc). (61) 

This supposition has been verified for g = 0.95 v^eV and various impurity concentrations. 
The values of remaining parameters correspond to (Lai_^.Ce^)Al2 alloy (Table The 
results, which are depicted in Fig. |9l confirm the hypothesis of the Jaccarino-Peter compen- 
sation effect in the considered theoretical model. According to Fig. [H the values of Tci{%) 
increases with % and after reaching a maximum at 7/ = "Hk, decreases and finally falls to 
zero. 

The value of an exchange field "Hj is expected to increase with increasing impurity con- 
centration, since the number of the magnetic moments of impurities, which are antiferromag- 
netically coupled to conduction fermions also increases. As a result the value of an external 
magnetic field l-L required to compensat the perturbative effect of "Hj increases with c. 

The Tc{H) graph depicted in Fig. [9^ resembles the dependence of superconducting tran- 
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H [T] H m 



FIG. 9. Tc{'H) graphs of (Lai_^Ce^)Al2 for g = 0.95 veV and under varying impurity concentra- 
tion: (a) X = 0.0010, (b) X = 0.0019, (c) x = 0.0028, (d) x = 0.0057. The values of the parameters 
5, go, g' , GqPf and M are collected in Table Hill 

sition temperature on dopant concentration p, exhibited by high-temperature cuprates, 
e.g. La2_:,.Sr^.Cu04^^'", YBa2Cu30/^ Bi2_^Pb^Sr2Ca2Cu30lo^^ Bi2Sr2_^.La^Cu06^^ and 
iron-pnictides^^. Tc{p) initially increases almost linearly in p and after reaching a maximum 
at optimal doping level popt, decreases and finally falls to zero. 

Analogous behavior of the superconducting transition temperature is observed, e.g. in 
CeRhln^^, CeCoIn^^, when the superconductor is under pressure. 

B. First order phase transitions 

In the case, of first order phase transitions, the assumption that the gap parameter A 
vanishes at the transition temperature does not hold. According to the results obtained in 

n 

Ref. Ill, the superconducting transition temperature Tc possesses three solutions (Td > Tc2 > 
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Alloy 


X 


c[%] 


M 


5 [eV] 


GqPf 


9' 


90 


9 W^] 






0.00 







0.3030 






0.00 






0.55 


2 




0.2890 






0.080 


LaCe 





1.00 


4 


0.01 


0.2800 


10/7 


0.006 


0.095 








o 










n inn 






2.00 


6 




0.2500 






0.105 




0.0010 




1 




0.2610 






0.10 




0.0019 




4 




0.2515 






0.189 


(Lai_a;Ce2,.)Al2 








0.01 




10/7 


0.006 






0.0028 




7 




0.2435 






0.19 




0.0057 




8 




0.2250 






0.23 






0.0 







0.2200 






0.00 


ThGd 




0.1 


3 


0.01 


0.2010 


2 


0.004 


0.11 






0.2 


5 




0.1915 






0.18 



TABLE III. The parameter values. 



Tcs) for certain values of g and c. These solutions can be determined numerically from the 
following equations: 

T,i : Fp{U) - Fscin) = 0, (62) 

Tc2: Fscin) -F^{'H) = 0, <^ = D,F, (63) 

Tes: Fpin) - FniU) = 0. (64) 

The existence of Tcs depends on the type of phase transition occuring at Tc2. If the system 
undergoes a phase transition to ferromagnetic phase at Tc2, then Tcs > for certain values 
of g. If Tc2 = TscD, then = and the system does not reenter the superconducting 
phase {SC or D). 

The solution of Eqs. fl62|) - flM|) for the parameter values corresponding to (Lai_a;Ce2^)Al2i 
are depicted in Fig. [101 For He (OT, 3.55T), three phase transitions P — )■ SC F ^ D 
with decreasing temperature are present. These phase transitions can be interpreted as 
NonSC — 7- SC —7- NonSC — > SC transitions, i.e. the Jaccarino- Peter compensation effect, 
which has been experimentally observed in a number of superconducting magnetic alloys, 
e.g.: Sn^.Eui.2-:.Mo6Ss^, Euo.TsSno.ssMogST.sSeo.g^^, CeCoIn^, URhGe i UCoGei^. 
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H[T] 



FIG. 10. Tc{n) graphs of (Lai_^Ce^)Al2 for g = 0.5 VeV and x = 0.0010. The values of the 
parameters 5, go, g' , GqPf and M are collected in Table Hill For V. G (0T,3.55T) the system 
undergoes three phase transitions with decreasing temperature, P — >• SC F ^ D, showing the 
Jaccarino-Peter compensation effect. 

VII. CONCLUDING REMARKS 

Our recent work on the phase diagrams of a BCS superconductor perturbed by a reduced 
s-d interaction^ was continued to examine the critical magnetic field and the effect of an 
external magnetic field "H on the superconducting transition temperature Tc of such system. 
The superconductivity enhancement, revealed by an increase of the critical magnetic field 
with decreasing temperature and increasing impurity concentration, has been found. 

Good quantitative agreement of the resulting theoretical expressions for the critical mag- 
netic field with experimental data was demonstrated for LaCe, ThGd and SmRh4B4. 

The numerical analysis of Tc{'H) showed that the perturbative effect of magnetic im- 
purities can be compensated by an external magnetic field, providing the evidence of the 
Jaccarino-Peter effect, which has been experimentally observed in a number of supercon- 
ducting magnetic alloys, e.g.: Sn2;Eui.2-a;Mo6Sg^, Euo.75Sno.25Mo6S7.2Seo.g^. The enhance- 
ment of superconductivity, displayed by an increase of the the upper critical magnetic field 
with decreasing temperature due to the Jaccarino-Peter effect was discovered in SmRh4B4— , 
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CeCoIngi^, URhGe and UCoGei^. 

The enhancement of the critical magnetic field in magnetic superconductors, which is 
due to the interplay of the superconductivity and magnetic order, opens the door to their pos- 
sible industrial applications. Since the upper critical magnetic field at absolute temperature 
of PbMoeSg magnetic superconductor, i.e. i/c(0) ~ 60T^ doubly exceeds the corresponding 
value of the high-performance, low-temperature superconducting material NbaSn, which is 
currently widely-exploited in power applications, e.g. in accelerators available at Fermilab, 
Brookhaven, DESY and CERN^. 
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